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Abstract. The Lefschetz-thimble approach to path integrals is applied to a one-site 
model of electrons, i.e., the one-site Hubbard model. Since the one-site Hubbard model 
shows a non-analytic behavior at the zero temperature and its path integral expression 
has the sign problem, this toy model is a good testing ground for an idea or a technique 
to attack the sign problem. Semiclassical analysis using complex saddle points unveils 
the significance of interference among multiple Lefschetz thimbles to reproduce the 
non-analytic behavior by using the path integral. If the number of Lefschetz thimbles 
is insufficient, we found not only large discrepancies from the exact result, but also 
thermodynamic instabilities. Analyzing such singular behaviors semiclassically, we 
propose a criterion to identify the necessary number of Lefschetz thimbles. We argue 
that this interference of multiple saddle points is a key issue to understand the sign 
problem of the finite-density quantum chromodynamics. 


PACS numbers: 03.65.Sq, 71.10.Fd, 71.15.-m 

1. Introduction 

Understanding strongly-correlated quantum many-body systems has been an ultimate 
goal in contemporary physics. Numerical simulation formulated on the discretized 
spacetime, especially lattice quantum Monte Carlo method, is a powerful ab initio 
tool for that purpose. Exact diagonalization of a Hamiltonian provides us complete 
information on the system; however, it usually requires the huge computational cost 
and is limited to small systems. Monte Carlo method, on the other hand, is based 
on the importance sampling in the phase space of the system, whose computational 
cost scales algebraically with the system size. Many thermodynamic quantities can be 
computed for various systems using this method, such as hnite-temperature quantum 
chromodynamics (QCD) in hadron physics [1], and liquid helium [2], ultracold atomic 
gases [3], Bose-Fermi mixtures [4] in condensed matter physics. 
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In many quantnm systems of great interest, however, Monte Carlo simulation suffers 
from the so-called sign problem [5, 6]. The path-integral weight e“‘^ becomes negative 
or even complex, and thus the importance sampling cannot be applied. Since the 
cancellation between large positive and negative noises gives a physically meaningful 
signal, the computational time to reduce statistical errors grows exponentially with the 
system size [7]. In condensed matter physics, famous examples of the sign problem 
are strongly-correlated electrons, such as the Hubbard model away from half-hlling [5], 
frustrated spin systems, such as the XY model on the Kagome lattice, and so on. In 
hadron physics, QCD at hnite quark densities attracts much attention for understanding 
the interior structure of neutron stars [8-10], but we have no ab initio simulation due 
to the sign problem [11]. 

The sign problem of the hnite-density QCD is known to become too severe when the 
quark chemical potential exceeds half of the pion mass [12-17]. At low temperatures, 
there is no phase transition until the quark chemical potential reaches about one third 
of the nucleon mass; however, the Monte Carlo simulation of two flavor QCD shows a 
phase transition at half of the pion mass [12, 13]. No one can describe the true onset 
at one third of the nucleon mass by using the path-integral formulation of lattice QCD 
[14, 15]. This is one of the biggest issue in hnite-density QCD. In this paper, we address 
this problem, which is called the “Silver Blaze problem” [14], by studying a toy model. 

There have been many attempts to attack the sign problem. The idea of 
complexihcation of the integration variables is recently developing, and the complex 
Langevin method [18-21] and the Lefschetz-thimble approach [22-25] to the path 
integral rely on this idea. The Picard-Lefschetz theory gives an extension of the steepest 
descent method to multiple oscillatory integrals, and Lefschetz thimbles correspond 
to steepest descent paths in the complex contour integration [22]. After hnding 
its usefulness in studying analytic properties of the Chern-Simons theory with the 
semiclassical analysis [23-25], it leads to diverse applications to real-time quantum 
phenomena [26, 27], resurgence trans-series methods of perturbation theories [28-33], 
and the sign problem in lattice simulations [34-42]. The Lefschetz-thimble approach 
is based on rigorous mathematics; however, its practical properties have not been fully 
understood yet. For recent developments, see Refs. [43-45]. 

In this paper, we apply the Lefschetz-thimble approach to the one-site model of 
electrons. This toy model can be regarded as an extreme limit of strong couplings 
because it can be obtained by neglecting hopping terms in the Hubbard model. The 
Hamiltonian of the one-site model can be easily diagonalized and thus we can calculate 
any expectation value exactly. However, since this model has the severe sign problem in 
its path-integral expression, it is hard to calculate expectation values by the conventional 
Monte Carlo method. This toy model provides us a good playground to study theoretical 
structures of the Lefschetz-thimble approach. In a previous study [34] , two-dimensional 
Hubbard model is studied using the Lefschetz-thimble Monte Carlo method, but one 
must use the so-called “one-thimble approximation” due to the current limitation of 
the numerical algorithm. Since our system is exactly solvable, we can understand 
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the complete structure of Lefschetz thimbles, and suggest appropriate approximation 
schemes. 

Based on the semiclassical analysis, we hnd that interference of complex phases 
among multiple saddle points is important to reproduce the step-function behavior of 
the density. The fermion spectrum at a complex saddle point resembles the quark 
spectrum of phase quenched hnite-density QCD, which would suggest the interference 
of multiple saddle points might also occur at hnite-density QCD. We will discuss this 
point in detail later. This study will help us to understand the sign problem in QCD. 

The outline of this paper is as follows. In Sec. 2, the one-site Hubbard model 
is explained, and its physical properties are calculated by exactly diagonalizing the 
Hamiltonian. The path integral expression of the toy model is introduced in detail, 
and the sign problem appearing in its path integral formulation is briehy reviewed. 
In Sec. 3.1, the Lefschetz-thimble method is introduced, and is applied to the one-site 
Hubbard model. Complex saddle points play a pivotal role in this approach, and thus we 
study their property in Sec. 3.2. We discuss structures of Lee-Yang zeros and the fermion 
spectrum at complex saddle points in Sec. 3.3. In Sec. 3.4, numerical results on the 
fermion number density are compared with the exact result using several approximations 
based on the Lefschetz-thimble approach. In Sec. 4, we discuss the resemblance of this 
model to the hnite-density QCD and address the Silver Blaze problem in QCD. We 
summarize our result in Sec. 5. 

2. One-site Hubbard model 

2.1. Physical properties of the model 

The Hubbard model [46-48] describes the lattice fermion system, whose dynamics is 
governed by the (second quantized) Hamiltonian, 

= -tYl + - h + H:i)- (1) 

* * 

The summation is taken over all the lattice sites i, and the notation {i,j) denotes that 
only the nearest neighbor hopping is considered, c^^and Co-,,, obeying = 

d^rdij, are creation and annihilation operators of fermions of the spin cr(='t', i) at the 
site i. is the number operator. The parameter t is called the hopping 

parameter, and we will consider only the case f = 0 in this section. The parameter 
U{> 0) describes the on-site repulsive interaction: The energy increases if two particles 
get together on the same site. The Hamiltonian is invariant under the global U{1) 
rotation i—>■ exp(±iQ;)c^^], and thus the total number of fermions is conserved. In the 
Hamiltonian (1), the chemical potential /i is introduced for this conserved quantity. 

In the strong coupling limit f = 0, one can solve the Hamiltonian (1) analytically. 
Let us see the result on the total number density {h) := {{h^ h^)) as a function of the 

chemical potential fi. Since each site is totally independent from others, one can study 
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H = — fan. 


( 2 ) 


In the following, we only consider the one-site case. Since the Hamiltonian (2) commntes 
with the number density operator h, we can take the number basis to hnd the ground 
state. Let us dehne the grand partition function by 


Z 


tr 


exp —jdH , 


(3) 


with the inverse temperature fd = 1/T. The number basis gives an explicit result: 


Z = l + 2e^^-t-e^(2^"^). 


(4) 


The number density is given as 




r\ = 1 7 = 

1 + 2e/5/^ + e/3(2M-f^) ■ 


(5) 


In the zero-temperature limit /? —)■ cxo, the number density takes 0, 1, and 2 for p/H < 0, 
0 < fi/U < 1, and pi/U > 1, respectively, and shows the step-function like behavior. In 
particular in the case pi/U = 1/2, we have (h) = 1 for any /3, which comes from the 
symmetry h i—>■ 2 — h (or, ho- e->• 1 — h^) in the Hamiltonian (2). This point is called 
half-hlling, at which the particle-hole symmetry exists. 


2.2. Path integral formulation and the sign problem 

We have explicitly seen that the one-site Hubbard model can be analytically solved by 
using the number eigenstates. However, if one changes the basis of the Hilbert space 
for taking trace, the sign problem emerges and thus the Hubbard model in the strong 
coupling provides us a good lesson. 

Let us hrst derive the path integral expression of the partition function (3). Using 
two-component complex Grassmannian variables the fermion coherent 

state is dehned by 

IV’) = exp(-#^)|0), (6) 

which satishes Co-lt/) = t/o-lV')- We take the convention in which all the Grassmannian 
variables anticommute with each other and with fermionic creation/annihilation 
operators. It is easy to check that 

(V’i|^2) = J drd^l^\^P)e-^*^{^P\ = 1, (7) 

and 

tr(0) = J d^jJ*d^jJe~'^*'^{—^jJ\0\^p). (8) 

Let In-i-, n^) be a fermionic Fock state, and then the inner products with coherent states 
are 

(0,0|^/) = 1, {l,0\'il)) = (0,1|^/) =-0^, (1,1|0) =-0t0t- (9) 
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Using these relations, one can compute the matrix element of the imaginary-time 
evolution operator as 


^exp. 


Ar, 


V'* (iji - 


( 10 ) 


Here we take an approximation for small Ar. In order to circumvent the quartic 
interaction, let us introduce an auxiliary real held (p via the Hubbard-Stratonovich 
transformation, that is. 


1 = \J^ j dv?exp-^(^(p-iU(^/>*e^"V)) • (H) 

Inserting Eq. (11) into Eq. (10), we obtain 

- \j ^ j dv? exp - + {‘ip- e^^'^(l + iArip)^'') ^ . (12) 

In order to take the naive continuum limit of the path integral expression, we need to 
exponentiate the auxiliary held in the fermionic bilinear term as 

e^^^(l + iAT(p) = exp (^Ar(/i + ip>) + ^{(p^)o + 

= exp (^Ar (^p+^ + icp^+ 0(Ar3/2)^ _ ^^ 3 ^ 

Here, {(p‘^)o means the expectation value of with the Gaussian weight in Eq. (12), 
which gives U/ At. The path integral expression of the partition function reads 



with the antiperiodic boundary condition ipN+i = —'f’l- Equation of motion in terms of 
(p shows 

(h) = ^ (15) 

The expectation value of (p is nothing but that of the total number density. This formula 
will be used later in order to compute the number density. Since Eq. (14) is quadratic 
in the fermionic held f) and "ip*, we can perform the Grassmannian integration explicitly. 
The spectrum of the fermion bilinear operator in the exponential in Eq. (14) is given by 

Xe{(p,p) = l-e(2^“d^dAfexp^^(iV3fc + /i + U/2), (16) 


and then, in the continuum limit, the fermion path integral (determinant) becomes 

=fl + exp/' dr {i(p{T) + p + U/2)\ . (17) 
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Since fermions have two spin degrees of freedom, the overall sqnare is taken. The fermion 
determinant (17) contains only the Matsnbara zero mode of (p. Thus the path integral 
of non-zero Matsnbara modes of (p gives a trivial Gaussian integration and does not 
depend on p. The path integral of our interest is now reduced to an integral of zero 
Matsnbara mode (pbg = Jq dr(p(r)//3, and we have 

Z = j d(pbg (l + (18) 

This integral can be performed analytically to hnd Eq. (4). Instead, we shall apply 
the Lefschetz-thimble method because Eq. (18) has the sign problem coming from the 
complex weight This simple model enables us to study the structure of the sign 

problem in terms of the Lefschetz-thimble approach. 

In order to explicitly show that the path integral (18) contains the sign problem, 
let us compare it with the phase quenched partition function 


Z' = 




2ttU 


d(pbg 


1 -I- e' 


,/3(i¥3bg+M+f^/2) 




This integration can be performed analytically as 


(19) 


Z' = 1 + 20®" + 


( 20 ) 


The positive sign in front of U is different from that in Eq. (4). If the last term is 
negligible, the sign problem is mild. This is the case for small chemical potentials 
Pl/U < —1/2 and low temperatures fdU 3> 1; the ratio Z/Z' is almost one because the 
dominant contribution comes from the unoccupied state. Therefore the sign problem 
is sufficiently weak. In the later section, we can interpret this semiclassically because 
the important conhgurations accumulate around the origin ipbg — 0. The sign problem 
becomes quite severe for p/U > —1/2 and low temperatures jSU S> 1. For example, at 
ja/U = 0, the ratio becomes 

Z/Z'~ 3e-^^ < 1. (21) 


This severe sign problem comes from the last term in Eq. (20), which causes the wrong 
onset of the density at p/U = —1/2. A similar problem happens in QCD at hnite 
density when the quark chemical potential exceeds the half of the pion mass [12, 13, 49]. 
For large chemical potentials, n/U > 3/2, we can switch the occupied and unoccupied 
fermions by changing the integration variable from ipbg to —p>hg + 2it/, so that the 
sign problem is weakened. There is another remarkable condition in which the sign 
problem disappears again by a simple shift of the integration variable: At the half- 
hlling p/U = 1/2, a new integration variable (pbg = -1- iU rewrites the original 

integral (18) into 



The integrand of Eq. (22) is nonnegative. This is because the particle-hole symmetry at 
the half-hlling becomes manifest in this new integration variable. As a conclusion of this 
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section, the sign problem can be circumvented in these three special cases, /i/f/ < —1/2, 
/i/T/ > 3/2 and fi = 1/2, even using the path integral expression; however, this does 
not happen in general cases. 


3. Lefschetz-thimble approach 


3.1. Lefschetz-thimble method 


Let us start with a multiple integration that gives the partition function 

Z= j (23) 

Jm" 

where S{x) is a complex action functional of the real held x = ... ,a;”) G In 

order to circumvent the oscillatory integral (23), we complexify the integration variables 
x^ H->■ z^ = x^ + iy^l. Integration is performed on steepest descent paths, called 
Lefschetz thimbles, instead of directly evaluating Eq. (23). Each Lefschetz thimble 
is an n-dimensional space spanned around a saddle point z^j in {a E T,). Let us 
consider Morse’s how equation for complexihed variables z by introducing a hctitious 
time t [22-24]: 


dz^t) ( dS{z{t)) \ 

dt V dz'^ J ’ 


Along this how, the integrand becomes smaller since 


—Re S{z{t)) = 


dS{z{t)) 


dz^ 


> 0 , 


(24) 


(25) 


and thus exp—S'(z(t)) —?• 0 as t —)■ oo. This how equation has an important conserved 
quantity: The imaginary part of S obeys 

^lmS{z{t)) =0. (26) 

Therefore, along each how line, the integrand has no oscillatory phase, which is a 
remarkable property to circumvent the oscillatory integral. Now, let us dehne a new 
integration cycle called Lefschetz thimble as follows: The Lefschetz thimble TTt is 
identihed as the set of points reached by some hows emanating from z^'. 


Ja — {^(0) I z{-co) ^ z^j. 


(27) 


Complex contour integral over jLa dehnitely converges since Re S diverges at boundaries 
of jTa-- On each integral, we have no sign problem because ImS* is constant. The 
partition function can now be computed by the sum of the nicely converging integrals 
as 

Z = d”ze-^(0. (28) 

o-es 


t There is another approach to the sign problem also based on the idea of complexihcation, the complex 
Langevin method [20, 21]. See Refs. [50-54] for recent developments. 
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The coefficient is given by the intersection number between and JC^- The dual 
thimble is defined as the set of the points reached by flows getting sucked into z^\ 

ICa = {2;(0) I z(+oo) = Za}. (29) 

There are two possible origins where the sign problem reappears in the Lefschetz- 
thimble method. One is a complex phase coming from the Jacobian of the integration 
measure d”z on Lefschetz thimbles, which is often called the residual sign problem [37- 
40]. Geometrically, the residual sign problem is mild if the Lefschetz thimble is almost 
straight. Another one comes from the possible cancellation in summing up Lefschetz 
thimbles in Eq. (28), because the complex phase ImS* can be different among several 
Lefschetz thimbles. 


3.2. Semiclassical analysis 


Let us apply the Picard-Lefschetz theory to the path integral (18). The effective action 
of this system is 

z^ - 2 \yi + exp (3 
The effective action satisfies the reality condition 

^ = ^(-^), (31) 



U 

iz + /i + - 


(30) 


and then the Lefschetz-thimble decomposition manifestly respects this symmetry so as 
to ensure the reality of physical observables [45, 55-57]. 

The logarithmic function has branch singularities at fermionic Matsubara modes 
z = i{fi U/2) + {2i 1 )'k/(3 for £ e Z. The integrand converges to 0, i.e., S —)■ J-oo, 
at these points. The flow equation of Eq. (30) reads 


dz _ /3_ 2i/?exp/3 + +f) 

dt 7/^'^i + exp/3(-i^ + /x + f)- 

Let us find the set of saddle points, which is an important step not only for the 
saddle-point approximation but also for the Lefschetz-thimble method. The saddle- 
point condition of the effective action (30) is 


2U 

1 -F exp -{3 {iZa + pi+^) 


(33) 


In order to simplify our analysis, we consider a limiting case where T <^U, \fi\. We can 
approximately obtain the saddle points as 

= i +T ^27rm + ihi + 0{T‘^) (34) 

for m G Z by assuming that the second term is much smaller than the first one. The 
reality condition (31) says that z^rn and Zm form a pair. If /i > 3[//2 or fi < —U/2, 
there is another solution 




2iU + o{T) for /i/G > |, 
0-|-o(T) for fi/U <—^. 


( 35 ) 
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Figure 1 . Behaviors of Morse’s downward flow equation for /3U = 30, t/ = 1, and 
n/U = 2 (a) [p/U = 0 (b)]. Star-shape black points show singular points of logarithm, 
and red blobs show complex saddle points z^- 


Behaviors of the downward flow equation (32) are shown in Fig. 1 with 17 = 1 
and (3U = 30. Let us discuss the case when the chemical potential is sufficiently large 
(/i/17 = 2). In this case, Morse’s flow is shown in Fig. 1 (a). Only one Lefschetz thimble 
77, associated with contributes in the Lefschetz-thimble decomposition. Other dual 
thimbles shown with dashed green lines do not intersect with the original integration 
cycle M as in Fig. 1 (a). Thus intersection numbers in Eq. (28) vanish and we have 


Z = [ dze-^^^\ (36) 

Jj, 

Within the saddle-point approximation using the complex saddle point z^, the number 
density is given as 


— 

n ~- 

U 


2 , 


(37) 


and thus the saturation of fermions is well described. The exactly similar thing happens 
also for /i/17 < —1/2, and then the saddle-point approximation gives n ~ —iz^/U = 0. 
This analysis explicitly shows why the sign problem is signihcantly weakened with an 
appropriate shift of integration variables for /i/17 < —1/2 and /i/17 > 3/2, as we have 
discussed in Sec. 2.2. 

In the following, we concentrate on the case where the sign problem is severe, 
—1/2 < /i/17 < 3/2, at low temperatures fdU > 10. In this parameter region, the 
typical behavior of the flow is shown in Fig. 1 (b). All dual thimbles /Cm intersect with 
the original integration cycle M. This shows that all the saddle points Zm contribute to 
the partition function, so that the interference among them may not be negligible. This 
interference requires a careful treatment of the semiclassical analysis in order to solve 
the sign problem. The big difference between Figs. 1 (a) and (b) can be explained by the 
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Stokes phenomena [23], which occur around n/U ~ 3/2 and also around ja/U ~ —1/2. 
However, this is irrelevant to the non-analytic behavior of the number density, and then 
we stick to the analysis for —1/2 < pi/U < 3/2. 

Let us denote the classical action at the saddle point as Sa{-= S{za)). 
Substituting the approximate expression (34) into Eq. (30), we have 

2 




(38) 


Re(^^- 

ImS'rn ^ 


2xm(/ + t 


(39) 

(40) 


For Eqs. (39) and (40), we have written down the m-dependent leading terms in terms 
of large f3U. Equation (39) shows that subdominant thimbles Jm can be comparable 
with the dominant one j7o for fdll S> 1 so long as |n7.|(7^ 0) is not so large. According to 
Eq. (40), these different thimbles have different complex phases, and thus a contribution 
from one Lefschetz thimble is generally canceled by other ones. Exactly at the half¬ 
filling, pi/U = 1/2, the complex phase is always an integer multiple of 27ri, and thus 
such a cancellation is absent. This gives an interpretation on the reason why the sign 
problem disappears at the half-hlling from the viewpoint of Lefschetz-thimble approach. 

Let us compute the partition function Zc\ only using these semiclassical information: 


Zc\ '■ = 



m=—oo 


(41) 


In this formula, Gaussian fluctuation around Zm is neglected because it only gives an 
unimportant overall factor at low temperatures. Indeed, within our approximation, 
diS{z^) is independent of m, and the result of the Gaussian integration slip by the 
summation. This summation (41) can be computed analytically by using the elliptic 
theta function as 

Zci ~ e-^° fl + 2 ^ cos 27rm 

\ m=l ^ ' 

= + (42) 

Here we used an approximate expression of the saddle points, Eqs. (39) and (40), to 
obtain Eq. (42). Using this result, in the limit (3 —>■ oo, the number density is given as 

(1 < fi/U < 3/2), 

(0 < fi/U < 1), (43) 

(-1/2 < /i/U < 0). 

This good agreement might be accidental to some extent, since we have neglected higher 
order contributions in the semiclassical analysis. This result still indicates the usefulness 
of the semiclassical approximation even when the sign problem is severe. 


1 5 

nd ■= InZci -> 

p op, 
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Let us analyze the non-analytic behavior (43) from the viewpoint of Lee-Yang zeros [58- 
66] in the semiclassical expression (42). Using the inhnite-product expression of the 
elliptic theta function, 

OO 

Hy ?) = n {} . (44) 

e=i 

We hnd that the zeros of the partition function in the complex pi plane satisfy 

cos (27r^) = cosh (^^(2^ - 1)^ (45) 

for some positive integer ^. Since Eq. (42) is originally derived as an approximate 
expression for —1/2 < pi/U < 3/2, we restrict the real part of the complex chemical 
potential to —U/2 < Re/i < 317/2. Then, the zeros of the semiclassical partition 
function are 


pi = ±i7rT(2£ — 1 ), pi = U ± mT{2i — 1 ). 


( 46 ) 


In the limit T —>■ 0, these zeros form straight lines crossing to the real axis at /i = 0, U. 
This behavior of zeros is a signal of the hrst-order phase transition in the context of the 
Lee-Yang zero analysis of a therodynamic system, and indeed we hnd the non-analytic 
jumps (43) of n at /r = 0,17 in the limit of T = 0 in the one-site Hubbard model. 

Let us apply the same analysis to the phase quenched partition function. We dehne 
the phase quenched semiclassical partition function by 



= ( 0 , 6 -"'’"’^') 


( 47 ) 


Note that this expression does not correspond to the semiclassical limit of Eq. (19) 
because the saddle points at —1/2 < pi/U < 3/2 do not he on the original integration 
cycle M. In the phase quenched approximation, the interference of complex phases 
g-iirnSm disappears, and then the /i dependence is given only by e“‘^°. Lee-Yang zeros 
do not exist if —U/2 Re/i 317/2. Therefore, the only possibility to explain 
the non-analytic behavior in this phase quenched approximation originates from the 
Stokes jumps around pi/U ~ —1/2, 3/2, which naturally explains the non-analytic but 
continuous change of n (see Ref. [44]). 

In order to investigate this failure of the phase quenched approximation more deeply, 
we analyze the fermion spectrum (16) at complex saddle points. For simplicity, let us 
take the continuum limit at hrst with an appropriate normalization, 

:= lim TNXe{ip,fx) = - |(2£ - l)7riT + iip + fx + ^\ . (48) 

N^cx) I / J 

The spectral representation of the number density at a background held ip is given by 

= -2T ^ ) = e-4(iv^+/4+c//2) + R 

i=—00 ^ \ f ! I / 


( 49 ) 
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If the sign problem is mild enough, p, < —U/2, then the eigenvalue at the saddle point 
2 ;* = 0 is 




- l)7riT + /i + 



(50) 


Since ReA^™(z*,p) cannot be zero, the chemical potential dependence of the number 
density is exponentially suppressed as, in the limit T —)■ 0, 


n{z^,p) ~ ^ 


~ 0 . 


(51) 


Q—y{fi+u/2) _j_ ^ 

However, for p > —U /2, Re A“° can be zero in the vicinity of saddle points, and indeed 
one hnds 




^(2(£ + m) - l)7riT-Tln 


^ + p ) 


(52) 


Within this approximation, the number density becomes 

+1 - (7 + )■ (53) 

The near-zero fermionic mode provokes the fake onset of the number density at 
p = —1//2 again. From the viewpoint of the fermion spectrum at complex saddle 
points, nothing special happens in the vicinity of the correct transition points p = 0,U. 
In the language of Lefschetz thimbles, this implies that their topological structure does 
not change at all around p = 0,U, as we can see in Fig. 1 (b). 


3.4- Numerical results 


Let us consider one-, three-, and hve-thimble approximations in order to analyze the 
importance of the interference among multiple Lefschetz thimbles. We consider the 
case where the sign problem is severe, i.e., —1/2 < p/U < 3/2 and [3U > 10. The 
(2m -|- l)-thimble approximation takes into account only Jq, J7±i,..., J±m in Fig. 1 (b). 
The partition function (18) in this approximation reads 




P m p 

/ dz ^ 2Re / d^ e'^F). 


( 54 ) 


Each contribution of the summation in the second term comes from a pair Nk and 571^, 
and becomes manifestly real because of the reflection symmetry ^ 1 —)■ —Since the 
Lefschetz thimble is homologically equivalent to the interval i{p U/2) -|- [{2k — 
l)7rT, {2k -|- l)vrT], the above integration (54) is reduced to 

/•+(2m+l)7rT 


^(2m+l) _ ^-S(x+i{fj.+U/2)) ^ 


( 55 ) 


— (2m+l)7rT 


with X = z — i{p U/2). The result on the number density at (3U = 30 is shown in 
Fig. 2, in which one-, three-, and hve-thimble approximations are compared with the 
exact computation. 





Lefschetz-thimble analysis of the sign problem in one-site fermion model 


13 



filU 


Figure 2. Behaviors of the number density n = —i{(p)/U as a function of g/U at 
j3U = 30. The solid black line shows the exact solution. Other lines, solid red one, 
dashed green one, and dotted blue one, show the result of one-, three-, and five-thimble 
approximations, which integrate over JTq, Uq U I7±i, and Jq U U J± 2 , respectively. 


The one-thimble approximation, which is shown with the solid red line in Fig. 2, is 
an approximation to integrate only over the Lefschetz thimble Jq. It almost gives the 
mean-held result, that is. 


h , 1 
riMF = y + 2- 


(56) 


for —1/2 < pi/U < 3/2. Figure 2 shows that the one-thimble approximation is 
not sufficient to describe the plateaus of the number density in the region —1/2 < 
pi/U < 3/2. One can easily check that Eq. (56) is also obtained from the phase 
quenched approximation (53) after shifting the integration variable in Eq. (18) as 
<Pbg </9bg -b i(/i -\-U/2). 

The result of the three-thimble approximation provides us a useful lesson for an 
application of the Lefschetz-thimble approach to the sign problem. The number density 
diverges at several chemical potentials around a rapid crossover, although the result is 
improved around each plateau. Let us analyze this divergence by using the semiclassical 
analysis. For that purpose, we introduce the semiclassical partition function with the 
(2m -b l)-thimble approximation by 


/^(2/Tl-t-l) 

^cl 



e 


-Sk 


k=—m 


(57) 
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In the three-thimble approximation, the semiclassical partition function approximately 
behaves as 


^_{3) ^ g-5o 


1 -h 2e cos27r f ^ -h - 


(58) 


Since fdU S> 1 and thus e ^ second term is not suppressed, compared with 

the hrst term. This approximate partition function vanishes at some chemical potentials 

/o\ 

around pi/U = 0 and 1 vanishes at pi/U ~ ±0.04, which is roughly consistent with 
Fig. 2). This example clearly demonstrates that picking up a part of the Lefschetz- 
thimble decomposition can violate the physical requirement, such as the positivity 
condition on the thermodynamic quantities. In this instance, the incompressibility 
{dnldp)lr? should not be negative. 

According to Fig. 2, we can conclude that hve Lefschetz thimbles JTqUJ 7±i U57±2 are 
necessary in order to explain rapid jumps of the number density in terms of the chemical 
potential at a low temperature ftU = 30. In fact, the hve-thimble approximation (dotted 
blue line) well reproduces the exact computation (solid black line) as seen in Fig. 2. How 
many Lefschetz thimbles are necessary in general at a given lower temperature jSl We 
propose the criterion to neglect the Lefschetz thimbles ffm for large |m|: 


7(2m+l) 


''cl 


-z. 


(2m—1) 


cl 


7(2m+l) 

^cl 


< 1 . 


(59) 


Assuming that the number of Lefschetz thimbles, (2m ± 1), is sufficiently large, one can 
replace in the denominator by the semiclassical partition function Zc\. Since 


“^cl 

I fy{27Tl —l) 

l^cl “ ^cl 


< 21 exp —S„ 
21 exp —S„ 


I, the contribution from can be negligible if 


< 


^cl 


(60) 


Here £ is a controlling parameter of an error, and e 1. Using the approximate results 
(39) and (42), we can solve this inequality with respect to m: 


m 


> 


/9U 8 

-- In -o's 

27r2 2 


TT 


U 1 
— ± - 
U 2 


^-27t^II3U 


(61) 


This criterion gives different results depending on /r, and thus let us first derive the 
strongest restriction. If the sign problem is severe, i.e., /i = 0, the elliptic theta function 
is exponentially small with respect to (3U: 


93 (^,e 


-2-k'^/PU 


2/Sf/ 


/3C//8 


TT 


(62) 


Therefore, for reasonable e such as £ ~ 0.1, the effect of £ is negligible and the 
criterion (61) gives 

fdU 


l“l ^ 4. 


(63) 


in the limit of f3U ^ 1. It means that we need at least {2\(dU/A ti'] ± 1) thimbles in 
order to describe the rapid crossover of the number density in the case of the one-site 
Hubbard model. 
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Let us discuss this behavior from another point of view. Since j3 becomes large, 
the relevant integration region becomes smaller as |Re z\ < ^U/13. However, Fig. 1 
shows that the length of each Lefschetz thimble is proportional to 1/(3, and thus we 
need more Lefschetz thimbles in order to cover the relevant integration region. That 
number is clearly proportional to yj(3U, and is much smaller than Eq. (63). Because of 
the interference of complex phases among Lefschetz thimbles, the semiclassical partition 
function becomes much smaller than that of the naive expectation if pi/U ~ 0, 1. This 
indicates that it is much easier to reproduce the plateau in the vicinity of the half-filling 
than that of the rapid crossover. The necessary number of Lefschetz thimbles for each 
phenomenon is also largely different. In fact, for pi/U = 1/2, one can show that 


.-So 




and, roughly speaking, the criterion (61) gives 


m 


> 


Ey/W 


(64) 


(65) 


for Jm being negligible with keeping a good accuracy only around the half hlling. This 
is consistent with the result of the above heuristic argument. Now, we can understand 
that the large gap of these two estimates (63) and (65) comes from the sign problem 
in summing up Lefschetz thimbles, and the semiclassical analysis provides a reasonable 
rough indication. 


4. Speculation on the Silver Blaze problem in finite-density QCD 

We studied the one-site Hubbard model as a toy model of the sign problem, and the 
Silver Blaze problem. We elucidate that the interference of complex phases among 
multiple saddle points is important to analyze these problems. In this section, we 
discuss the analogy of the Silver Blaze problem in the one-site Hubbard model with 
that in finite-density QCD. 

Let us pay attention to the behavior of the baryon number density hb for finite- 
density QCD at T = 0. In general, the quark determinant, Det [ 74 (/^a + "^) ~ A^qk], 
will have the /iqk-dependence, which naively means that the baryon number density 
arises for any chemical potentials. However, since we know that the lightest baryon 
in the QCD spectrum is nucleon, which is a composite particle of three quarks, the 
baryon number density hb must be zero for any 3/iqk < — B ~ 923 MeV (mjv is 

the nucleon mass and B is the nuclear binding energy). The /iqk-dependence of the 
QCD partition function must vanish for /iqk < m^yf/S at T = 0, which empirical fact is 
highly nontrivial from the viewpoint of the path integral of lattice QCD as discussed 
above (Fig. 3 illustrates this situation schematically). Understanding this behavior of 
Hb must be a key first step to correctly perform the lattice QCD simulation at arbitrary 
temperatures and baryon chemical potentials. In Ref. [14], this is named the baryon 
Silver Blaze problem after the famous story of Sherlock Holmes. 
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Figure 3. Schematic illustration of the Silver Blaze problem in the finite-density 
QCD. The baryon number density jumps at /r ~ mjv/3 (solid blue line). Phase 
quenched theory shows the early onset of the baryon number density (dashed red 
line). 


So far, the baryon Silver Blaze problem has been well understood only for pqk < 
m7r/2 (m^r — 140 MeV is the pion mass) [14, 15, 65]. In general, the eigenvalues 
of 74 (.0A + m) — /iqk can be expressed as A(j,n)(y4,pqk) = £j{A) — pqk — + 

i(2?7, + l)7rT with —vrT < (fj < ttT [14], Therefore, Det [74(.0^ + m) — Pqk] = 
X(j,n)- The product over n can be taken in the exactly same way as we 
usually do for that over Matsubara frequencies, and we hnd that it is proportional to 
_j_ g-/3(ej-i</-:/-Mqk)/2j igPg} 'jpe fermion determinant becomes 

Det [ 74 A + m)- pqk] 

= A/"]^ _j_ g-/thi-i</'j-Mqk)/2j 

j 

= J\f n(i+ g-/3(£j-Alqk)gWj j 


Vff 


Afe'+i'/'j) 


(1 + 


^ /^qk '^4 


7(1 + 


i“/5(si+A‘qk+i7j) 


). 


( 66 ) 


£,>0 


where M is the normalization factor (see Ref. [15] for its more rigorous derivation using 
the zeta-function regularization). Its /iqk-dependence at sufficiently low temperatures 
and at Pqk > 0 becomes crystal clear according to the following expression [14, 15]: 

Det [ 74 A + m)- /iqk] 


Det [74(.0 a + rn)] 


(24.)<C/^qk 


(67) 


The quark determinant becomes independent of /iqk at T = 0 if /iqk is smaller than the 
minimum of ej{A){> 0). For statistically signihcant gauge helds, it is conhrmed that 
min(ej(il)) = m^/2 [14, 15, 65], and thus = 0 for /iqk < This also means that 

the Lefschetz-thimble decomposition of lattice QCD at /iqk < ni7r/2 and at sufficiently 
low temperatures is identical to the original integration cycle up to an exponentially 
small correction because the quark determinant at /iqk = 0 is positive. 

In contrast, as computed in Eq. (17), the fermion determinant in the one-site 
Hubbard model is 


Det 




i(p(r) + h + Y 


= (1 + 




) 


( 68 ) 
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We can ask the same question: How can we show n = 0 for /i < 0 in the limit (3 = oo 
by using the path-integral expression when the fermion determinant depends on /x? 
This Silver Blaze problem for /i < —U/2 can be solved in the same manner as the 
baryon Silver Blaze problem at /Xqk < m7r/2. At /i < —U/2, the fermion determinant 
becomes independent of /i as T —)• 0. This also explains why the Lefschetz thimble 
jZi, at /i < —U/2 is almost identical to the original integration cycle M (The similar 
discussion shows the Lefschetz thimble at /r > 377/2 is almost identical to the line 
with Im(z) = 2). 

For —U/2 < /i < 0, the /x-dependence of the fermion determinant (68) becomes 
exponentially large, and thus the fact that n = 0 at T = 0 is still veiled. At /x ~ —U/2, 
zeros of the fermion determinant (star-shape black points in Fig. 1) come close to the 
real axis, which is almost identical to the Lefschetz thimble for /x < —U/2. In 
the case of the one-site Hubbard model, this triggers Stokes jumps and the original 
integration cycle is decomposed into multiple Lefschetz thimbles for /x > —U/2. One 
can no longer replace the expectation value of the number density by the number density 
at a complex saddle point. The significance of interference among multiple Lefschetz 
thimbles is identihed in order to explain not only the rapid jumps of n but also the /x- 
independence of rx (n = 0) for —U/2 < /x < 0 at the zero temperature. Furthermore, if 
the number of Lefschetz thimbles is insufficient, the thermodynamic stability is violated 
in the one-site Hubbard model. 

We speculate that this decomposition of Lefschetz thimbles and accompanying 
interference of complex phases also play an important role in finite-density QCD, in 
particular, in the Baryon Silver Blaze problem and in the sign problem beyond half of 
the pion mass. For /Xqk > mT,/2, the quark determinant becomes dependent on /Xqk and 
highly oscillatory on the statistically significant domain of real gauge fields. This means 
that zeros of the quark determinant come close to that domain. It would be natural to 
conclude that the statistically dominant region of the original integration cycle starts 
to be decomposed into multiple Lefschetz thimbles around /Xqk — mT,/2. The complex 
phase of the fermion determinant becomes ill-defined at its zeros. Thus, the constant- 
phase condition of the integrand on a Lefschetz thimble would be hard to be satisfied 
without such decomposition, when zeros of the fermion determinant are located in the 
vicinity of the Lefschetz thimble. Once the original path integral is decomposed into 
multiple thimbles, we might need to carefully take their interference into account to 
explain the /Xqk-independence of until /Xqk exceeds {mN — B)/3. 

Because of the complexity of QCD, we cannot show these statements rigorously 
so far. Future study of QCD-like models based on Lefschetz thimbles and justifying 
our speculation will be crucial to develop our understanding on the baryon Silver Blaze 
problem and on the sign problem of the finite-density QCD beyond half of the pion 


mass. 
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5. Summary 

One-site repulsive Hubbard model shows a non-analytic behavior by changing the 
chemical potential at the zero temperature (3 = oo. This model is obtained by taking 
the strong coupling limit, and thus our analysis of Lefschetz thimbles forms basics of 
studying the strongly-coupled repulsive Hubbard model using the path integral. Since 
the number eigenstate can solve this model exactly and its path integral expression has 
the sign problem, this is a good testing ground to analyze the structure of the sign 
problem in various solutions. Furthermore, the non-analytic behavior of the number 
density looks quite similar to the Silver Blaze problem in hnite-density QCD, and thus 
the understanding of this model from the functional integral itself is an important task. 

In this paper, the sign problem in the Lefschetz-thimble approach is carefully 
studied. Topological structures of the Lefschetz thimbles are analyzed based on the 
semiclassical analysis, and their relation to the fermion spectrum is investigated. If 
one picks up only the most relevant Lefschetz thimble, the mean-held approximation is 
almost recovered. The signihcance of interference among multiple Lefschetz thimbles is 
identihed in order to explain the non-analytic behavior. As we lower the temperature, 
there exist regions where the necessary number of Lefschetz thimbles increases linearly 
in jdU. On the other hand, in order to explain the Silver Blaze like behavior, that 
number only grows as ^/[3U. 

If the number of Lefschetz thimbles is insufficient, we found not only large 
discrepancies from the exact result, but also thermodynamic instabilities due to artifact 
of the approximation. This means that the Lefschetz-thimble decomposition does 
not manifestly ensure the thermodynamic stability, or the positivity condition on the 
partition function, although its reality is always true with a reasonable condition [45]. 
Developing our understanding on the stability from the viewpoint of Lefschetz thimbles 
is important. 

Since the one-site Hubbard model is quite a simple model, it is not clear whether 
interference of multiple Lefschetz thimbles is relevant for studying thermodynamic 
systems with complicated interactions. Nevertheless, we discuss the speculation on the 
baryon Silver Blaze problem of hnite-density QCD based on its mathematical similarity 
to the one-site Hubbard model. We conjecture that the statistically dominant gauge- 
held conhgurations are decomposed into multiple Lefschetz thimbles when the quark 
chemical potential exceeds half of the pion mass. To justify our conjecture on hnite- 
density QCD and QCD-like models is signihcantly important to solve the baryon Silver 
Blaze problem and to numerically simulate the hnite-density QCD beyond the half of 
the pion mass. 

In order to understand strongly-correlated electron systems with the sign problem, 
applying the Lefschetz-thimble method to the two-site Hubbard model is a nice 
straightforward extension of this study. This enables us to take into account the ehect of 
the hopping term. If the atomic potential is sufficiently strong, localization of fermions 
is a natural scenario and thus two-site models will become benchmark. We believe that 
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semiclassical analysis with complex saddle points provides ns a deeper nnderstanding of 
the sign problem, and developing that techniqne in general cases is an important fntnre 
stndy. 
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